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ABSTRACT 


The  reaction  between  water  and  sodium,  followed  by  reaction  of  the 
liberated  hydrogen  with  air,  in  a  closed  vessel  is  treated  as  a  thermodynamic 
problem  at  constant  volume  „  The  results  for  reaction  of  stoichiometrically 
equivalent  amoimts  of  sodi,um  and  water  for  varying  size  of  vessel  are  found 
to  follow  the  equations  P  =  3o44  +  45»64  R  “  7.97  R^  where  P  is  the  final 
pressure  in  pounds  per  square  inch  gauge,  and  R  is  the  ratio  of  sodium  (in 
hundreds  of  cubic  feet  measured  at  350*^F)  to  dry  air  (in  millions  of  cubic 
feet  at  77”F) .  This  equation  is  satisfactory  when  the  final  ten^jerature  is 
above  the  melting  point  of  sodium  hydroxide  (591°K). 
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STATIC  PEESSUEE  DUE  TO  SCDIDM-WATEE 
EEACTIONS  IN  CLOSED  VESSELS 

W„  So  Horton 


IHTRCDUCTION 

In  the  development  of  the  SIR,  the  bringing  together  within  the  same 
general  area  of  large  quantities  of  molten  sodium  (as  coolant)  and  water  (as 
shield)  has  necessitated  an  inquiry  into  the  hazard  associated  with  the  acci” 
dental  mi  yi  ng  of  these  two  substances „  It  has  been  stated  that  at  reasonable 
distances  from  the  source  of  explosion  ary  structure  that  will  lathstand  the 
static  equilibrium  pressure  will  also  withstand  the  dynamic,  transient  ex¬ 
plosion  pressrure  which  generally  is  of  the  order  of  twice  as  great  as  the 
static  equilibrium  pressure.  It  appeared  necessary  therefore  to  compute  the 
static  equilibrim  pressiire  from  reaction  of  sodium  and  water  in  a  closed 
vessel.  For  the  calculations  described  herein,  it  was  desired  to  compute  the 
pressure,  the  temperature  of  the  gas  phase,  and  the  amount  of  energy  released 
or  available  to  elevate  the  temperature  of  the  materials  above  25°C,  Other 
information,  generally  provided  during  the  course  of  the  calculations  but  not 
described  herein,  includes  the  water  content  of  the  gas  phase,  the  heat  of 
hydration  of  sodium  hydroxide,  and  the  minimum  size  of  the  vessel  which  would 
ensure  that  the  pressure  would  be  no  greater  than  a  prescribed  amount.  This 
report  gives  a  detailed  explanation  of  the  similar  calculations  made  by  the 
author  for  hazard  evaluation  of  the  SIR. 

The  numerical  design  data  necessary  for  the  calculation  consisted  of s 

Ic  The  weight  of  sodium  involved, 

2.  The  average  temperature  of  the  sodiumo 

3.  The  weight  of  water  involved. 

4o  The  average  temperature  of  the  water. 

5o  The  volume  of  the  containing  vessel. 

6.  The  composition  of  the  original  gas  in  the  containing  vessel. 

7.  The  average  temperature  of  the  original  gas  in  the  containing  vessel. 

The  assumptions  made  in  the  calculation  include  the  followings 

1.  The  evolved  or  available  heat  is  partitioned  among  all  products  of 
the  reaction  but  not  among  any  other  materials, 

2.  hydrogen,  produced  by  the  action  of  sodium  on  water,  reacts  with  any 
oxygen  available  in  the  containing  vessel  and  liberates  heat. 
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3.  AIT  gases  behave  like  perfect  gases  under  the  conditions  of  the 
reaction o 

4o  The  volume  of  substances  destroyed  in  the  reaction  is  insignificant 
compared  to  that  available  in  the  containing  vessel. 

Ebjperiments  have  shown  that  the  sodium-water  reaction  is  not  very  fast  as 
compared  to  a  high  ejcplosive^  and  that  a  significant  fraction  of  the  liberated 
heat  is  given  to  neighboring  materials.  It  was  felt  wise,  in  this  calculation, 
however,  to  make  Assumption  1  for  a  number  of  reasons.  It  is  almost  impossible 
to  forecast  accurately  what  fraction  of  the  heat  will  be  so  distributed.  The 
answer  will  depend  upon  where  in  the  system  the  nixing  occurs  and  upon  the  un¬ 
known  time  constant  of  the  reaction  as  it  may  occur.  By  making  the  assumption 
a  safety  factor  is  provided,  since  accounting  for  the  neglected  heat  loss  would 
reduce  the  gas-phase  temperature  and  consequently  also  reduce  the  pressure. 

The  general  method  of  attack  applies  simple  chemical  thermo<^mamics  to 
the  two  reactions,  sodim  plus  water  and  hydrogen  plus  ojjygen.  The  heat 
liberated  at  constant  volume  is  partitioned  among  the  products.  The  final 
temperature  and  the  pressure  of  the  gas  is  computed.  The  method  is  similar 
to  that  described  by  (Hasstone.  The  main  difference  between  that  de¬ 
scription  and  the  one  in  this  report  is  that  here  constant  voliame,  rather 
than  constant  pressure,  is  assumed.  If  the  computed  pressure  is  greater  than 
the  bursting  pressure,  the  vessel  will  rupture.  If  the  reverse  is  true,  the 
computed  pressure  will  be  the  time  pressure  at  equilibrium  under  the  assumptions 
given  above. 

It  must  be  recognized  that  the  details  of  calculation  depend  upon  the 
physical  states  of  the  water  and  sodium  hydroxide  at  equilibrium.  Some  combi¬ 
nations  are  listed, 

1.  Condensed  water  region 

2.  Vaporized  water  region 

3.  Molten  sodium  hydroxide  region 

4.  Vaporisied  sodium  hydroxide  region 

5.  !^drated  sodium  hydroxide  region. 

Within  a  region  the  pressure  (or  temperature)  as  a  function  of  sodium  weight 
per  vessel  volume  ratio  (or  sodium  per  water  ratio)  is  continuous.  At  the 
boundary  between  two  regions,  however,  there  is  a  point  of  discontinuity  caused 
by  differences  in  the  heat  capacity  of  the  two  phases  of  the  changed  substance 
(e.g.,  water  vapor  versus  liquid  water)  and  the  heat  of  transition. 


11 


THF.RMQnYNAMC  ANALYSIS 

The  procsess  involved  may  be  represented  by  the  equation  (the  thermo— 
cynamic  notation  of  Lewis  and  Randall^^^  is  used  in  this  report) ; 

Na(j^,  6A4°K)  +  X  H20(j^,  339°K)  +  ^  C^(g,  2980K)  +  ^  %(g,  298°K)  =  NaOH(^,  T2) 

,  \  /  \  (1) 

+  f  Z  _  1  )  H20(g,  T2)  +/  £  -  i  C^(g,  T2)  +  ]^(g,  T2). 

yx  2  j  ^  lx  4y  X 

Here  the  sodium-water  and  hydrogen-oxygen  reactions  have  been  combined.  The 
indicated  temperat\ires  of  reacting  substances  apply  to  the  SIR  under  operating 
conditions.  Purej  dry  air  at  room  temperatiore  has  been  used  for  the  original 
gas  in  the  containing  vessel.  The  latter,  of  course,  is  a  buildingo  The  re¬ 
action  would  involve  x  gram-moles  of  sodium  with  y  gram— moles  of  water  in  the 
presence  of  5z  gram-moles  of  air.  The  equation  has  been  divided  by  x  since 
the  resulting  temperature  and  pressure  depend  upon  the  relative  quantities  and 
this  division  reduces  the  size  of  figures  carried  along  in  the  calculations. 

For  most  of  the  results  used  in  the  final  phases  of  the  SIR  hazards  evalu¬ 
ation,  the  molten  sodium  hydroxide  region  was  involved.  This  made  it  convenient 
to  adopt  591°K,  the  melting  point,  as  a  base  temperature.  The  desired  calcu¬ 
lation  goes  through  the  following  steps; 

1.  Determine  -  AE  (the  heat  evolved  at  constant  volume)  for  the  reaction 
with  591°K  for  T2. 

2.  From  the  combined  heat  capacities  at  constant  volume  of  the  products 
(including  oxygen  and  nitrogen)  and  AE,  determine  AT ,  the  actual  amount  by 
which  T2  is  above  591°Ko 

3 .  Compute  T2  =  591  +  At , 

4.  Compute  the  pressure  from  the  perfect  gas  law. 

Each  of  these  steps  will  now  be  treated  in  detail.  In  order  to  fix  the 
principles  more  firmly,  a  specific  example  will  be  carried  through  siraul- 
taneous.ly. 

The  assumed  example  will  concern  the  resifLt  of  mixing  10760  pounds.of 
sodium  with  a  chemically  equivalent  amount  (Z  =  1)  of  water  in  a  2  x  10°  cubic 
foot  vessel  containing  air.  This  leads  to  the  following  pertinent  data. 

Sodium  =  2,12  x  10^  gram-atoms 
Original  oxygen  =  O.464  x  1.0°  gram-moles 
Original  nitrogen  =  1,86  x  10^  gram-moles 


12. 


^  =  1,00  (assumed) 


S,  ~  0 » X  10^  =  o  ]  Q 

X  2,12  X  105“  ^ 


=  4  X  2,19  =  8,76 


X 

Z 

1^ 

1  =  0.50 

X 

2/ 

z 

_!) 

=  1.94 

X 

4  1 

To  cany  out  Step  Ij.  +.he  ste^iidard  values  for  AH  of  formation  as  obtained 
from  reliable  soiirces^^j  are  used  to  compute  for  the  reaction  when  all 
substances  are  in  their  staiadard  states  a,t  one  atmosphere  pressure  and  25°C, 
This ^ yields  AH®  =  -67 ,84  kilocalories „  To  change  this  to  constant  voluine 
conditions  (AE)  rather  than  constant  pressure  (Z!1I)  the  equation 


AE  =  AH  ^  (AN)  itr  (2) 

is  used  where  AN  is  the  change  in  the  number  of  raoles  of  gases  for  the  reaction. 
When  all  the  substances  are  in  their  respective  standard  states,  only  the 
oxygen  and  nitrogen  are  gaseous.  Therefore 

AN  =  i, 

4 

-(AN)  RT  =  -  (^0,25)  X  1.987  x  298  =  +15  cal  =  0.15  kcal, 
and  AE  =  -67.84  +  0.15  =  “67,69  keal. 


To  change  the  substances  to  their  respective  reaction  temperatures  and  states, 
the  following  transformations  are  added. 


(1) 

Standard  reaction 

AE  =  -  67,69 

(2) 

Na(jJ,  644)  =  Na(.^,  371) 

“  1.90 

(3) 

Na(i,  371)  =  Na(s,  371) 

-  0,63 

(4) 

Na(s,  371)  =  Na(s,  298) 

~  0.52 

(5) 

H20(/,  339)  =  H20(j^,  298) 

-  0.72 

(6) 

Na0H(s,  298)  =  Na0H(s,  591) 

+  7.47 

(7) 

Na0H(s,  591)  =  Na0H(i,  591) 

+  lo6o 

13 


(8) 

1/2  H20(i,  298)  =  1/2  H20(/,3T3) 

+  0.68 

(9) 

1/2  I^0(/,  373)  =  1/2  H20(g,  373) 

+  4.49 

(10) 

1/2  H20(g,  373)  =  1/2  H20(g,  591) 

+  0.92 

(11) 

1o94  02(g,298)  =  1.94  02(g,591) 

+  3.05 

(12) 

8.76  N2(g,298)  =  8.76  N2(g,591) 

+  11^ 

AE  =  -  40.22  kcal 

All  energies  in  thivS  set  of  12  transformations  are  calciilated  from  heat  ca¬ 
pacities  at  constant  volume  or  heats  of  transition  at  constant  volume,  ftom 
this  step  -  AE  =  40,220  calories  per  gram-atom  of  sodium  is  obtained — this  is 
the  heat  available  to  heat  the  products  above  591® , 


In  order  to  perform  Step  2  conveniently,  mean  heat  capacities  are  used 
(see  Reference  1,  pp.  52-4)'=  Available  equations  (Reference  1,  p„  503)  for  o 
were  converted  to  the  desired  quantities.  These  are  listed  here.  ^ 

liltrogen  (g)s  =  4.906  +  0.625  x  10“^  T2  cal/deg/g-mole  (3) 

Cb^rgen  (g)  s  =  4a957  +  1.369  x  lO^^  T2  -  0.3077  x  10“^  T2^  (4) 

E2O  (g)3  =  5.981  +  1.205  X  10"^  T2  +  0.943  x  10“^  12^  (5) 

The  heat  capacities  listed  are  mean  values  between  591°K  and  T20 
For  liquid  sodium  hydroxide  the  equation 

NaOH  {l^)%  =  18.22  -  3c836  x  lO-"^  T2  +^4^1  (5) 

'^2 

was  derived  from  that  given  for  NaOH  (£.)  ly  Epstein  and  Weber.  ^4)  These  four 
relations  are  plotted  for  various  values  of  T20  It  is  now  possible  to  deter- 
mine  T2  by  successive  approximations.  The  total  heat  capacity  is  given  by 


e  =  'e^ 


(7) 


where  c^  is  the  mean  heat  capacity  of  the  1th  product ^  n;j^  is  the  coefficient 
of  that  product  in  the  balanced  chemical  equation  (l) .  The  difference  between 
Tg  and  591°K  is  given  by 


(8) 


and  T2  then  Iss 


T2  =  591  +  at 


(9) 


14 


Now  if  the  ej_’s  are  approximated,  an  Initial  approximation  for  T2  can  be 
calculated  from  the  combination  of  Equations  (7),  (8),  and  (9)  or 


T 


2 


=  591  + 


•-AE 

2%  Ci 


(10) 


This  gives  a  new  value  for  T2.  The  process  should  then  be  repeated  until  no 
change  occurs  witliin  the  desired  accuracy .  For  the  numerical  example  Included 
here  the  following  table  may  be  set  up  for  the  initial  approximation  of  T20 


Substance 

Ei 

ni  Ci 

NaOH(i?) 

1  18  ( 

cal/mole/ deg 

18  cal/deg 

H20(g) 

0,50  8 

4 

OgCg) 

1,94  6 

12 

»2(g) 

8.76  6 

52 

C  : 

=  2n^  Ci  = 

86 

Prom  this  we  obtain 

AT  =  4^0  =  ^58^ 

86 

and  T2  =  591  +  468  =  1059°K. 

From  the  graphs  of  the  new  estimates 

repeated  as  indicated  below. 

are  obtained  and  the  calculation  is 

Trial 

c,  cal/deg 

AT, 

T2,  °K 

1 

86 

468 

1059 

2 

82,58 

487 

1078 

3 

82,60 

487 

1078 

The  final  step,  computation  of  the  pressure,  requires  the  molar  volume  of 
the  gases ,  The  total  niimber  of  moles  of  gas  is  given  by 

N  =  S 

(11) 

where  S  is  the  number  of  gram-atoms  of 
summation  over  gaseous  products  only. 

sodi’jm  and  the  "prime"  indicates 

It  should  be  noted  that  the  air  in  the 

vessel  is  considered  to  be  a  product  since  it  has  been  factored  into  the  coef¬ 
ficients  of  Equation  (l) o  Thus  for  the  chosen  example; 

N  =  2ol2  X  105  (0,50  +  1.94  +  8.76)  =  2,374  x  10^  moles 
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The  molar  -yolume  is  gi,Yen  by? 


(12) 


where  V  is  the  volume  of  the  containing  vessel  In  liters »  For  the  example 
again, 

V  =  =  23  086  liters/mole 

2e374  X  10®  moles 


The  pressure  is  now  given  by  the  perfect  gas  l.aw 

p  =  15p_ 

V 


(13) 


where  R  is  the  mcxlar  gas  constant o  The  example  gives 


P 


(0._Q82Q6lilQ78l 

23«86 


3c7a7  atm„ 


The  remaining  figure  required  l.s  the  energy  av.ailable  to  heat  materials 
above  25'^C,  Tiiis  l,s  obtained  by  the  summation  of  reactions  on  page  12  using 
only  the  first  five  reactions.  This  gives  for’  one  gram~atom  of  sodim  a 
AE  of  -71o48  kcal„  Q,  the  total  heat  evolved  at  25° ?  then  is  obtained  from 
the  equation 

Q  =  S  (-AE  298) 

=  2„12  X  ICh  (7c.148  X  104) 

=  lo;52  X  10^’*^  cal.orles 


This  development  can  be  used  for  any  phase  region  uti.lizing  the  lower 
temperature  of  the  region  as  a,  basCc 

DEFELOIMEICT  OF  EQMTIOM  AxND  GEI^RAIiZED  GWES 

The  method,  althoiigh  simple,  requires  time  and  when  a  large  number  of 
random  points  ai'e  requested,  it  was  useful  to  develop  a  more  direct  calcu¬ 
lation  reqi-urlng  but  a  fe'w  minutes  „  After  a  number  of  such  points  were 
supplied  a  method  was  sought  for  representing  the  results  in  an  equationo 
It  appeared  at  first  that  a  complex  equation  would  be  required,  with  quantity 
of  sodium  and  volume  of  container  as  separate  independent  variables.  However, 
it  was  felt  that  by  using  'the  ratio  of  sodiaim  to  container  volume  as  a  single 
variable  a  satisfactory  equation  could  be  developed.  The  reasoning  follows 
the  steps  of  the  nuroBrlcal  analysis  just  given. 
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If  the  steps  for  the  calculation  of  -AE,  hereafter  to  be  called  Q,  are 
examined  it  is  seen  that  this  quantity  can  be  expressed  as 


Q  = 


E, 


(15) 


This  conresponds  to  dividing  the  energy  into  parts  due  to  all  other  substances 
(Reactions  1  through  lo)  and  those  due  to  the  air  of  the  container  (11  and  12) . 
The  ratio  of  water  to  sodium  is  considered  constant  and  will  be  discussed  later 
in  the  report.  The  equation  can  be  rearranged  as 

s  =  %  -  )  (16) 


where 

and 


Ee  =  %  +  i 

E'  =  |(4£b  +  Eo) 


Since  ^  is  the  ratio  of  air 

X 


to  sodim,  it  is  replaced  by  R  giving 


Q  =  -  E’R. 


(17) 


This  shows  that  the  energj^  is  a  linear  function  of  the  air-to-sodium  ratio. 


The  combined  mean  heat  capacity  of  the  products  may  be  xjritten  as 

C  =  f^  (T)  +  i  f2  (T)  + /z  _  1  (T)  +  4z  f,  (T)  (18) 

2  I  X  4-  /  X  ^ 

where  the  functions  of  T  are  given  for  sodium  hydroxide,  water,  oxygen,  and 
nitrogen  in  Equations  (3),  (4)>  (5)  and  (6),  They  can  be  rewritten  as 

G  =  f(T)  +  ^  g(T) 

X 


or 


C  =  f (T)  +  g(T)  »  R 


The  expression  for  AT  now  becomes 


AT 


T  -  Tr  = 


Ej^  -  E'R 
f(T)  +  g(T)R 


(19) 

(20) 


Since  the  functions  f(T)  and  g(T)  are  at  least  quadratics.  Equation  (20)  becomes 
a  cubic  if  a  solution  is  sought  for  T  in  terms  of  R,  It  appears  reasonable  to 
assume  that  an  explicit  expression  for  T  could  be  found  by  expanding  in  powers 
of  R,  Thus  Equation  (20)  is  replaced 


T  =  Aq  +  Bi  R  +  C^R^  + 


o  e  o 
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To  continue  the  development  the  molar  volume  is  needed  „  The  total  moles 
of  gas  is  given  by  {0c5  +  2.  „  1  or  (0o25  +  The  total  volume  en- 

X  4  X  _ 

closing  the  gases  is  V  or  equal  to  Vg^(5z)  where  is  the  molar  volimie  of  air 
at  normal  conditions.  The  molar  volume  of  the  product  gases  then  becomes 


F  -  5^¥a  R 

“  t0o25  +  “  (0.25  +  R) 


(21) 


The  computed  absolute  pressiure  according  to  the  perfect  gas  law  is; 


P  = 


(0.082) 


At  +  BtR  +  C^R^  + 


0»082  /  0.25 (At  +  B^R  +  GjR^  + 
Ta  V  R  ) 


ft  o  & 


) 


Thus  P  =  f(R)  (22) 

Since  it  is  now  shown  that  the  pressure  is  a  function  of  R,  the  air-to- 
sodium  ratio,  it  is  evidently  a  ratio  of  R*  =  R“^,  also.  The  latter  proves 
more  convenient  to  xise.  It  seems  reasonable  to  take  a  number  of  values 
already  computed  by  the  long  procedure  and  determine  an  equation  of  the  form 

P  —  a  +  bR*  +■  cR®  ^  +  ... 

to  fit  the  data.  The  remaining  problem  now  appeared  to  be  the  number  of  terms 
for  a  satisfactory  fit.  Actually^  the  initial  attempt,  which  was  for  an  approxi¬ 
mate  equation,  proved  to  give  the  desired  accuracy  when  checked  against  all 
previously  calculated  points.  Thi.s  was  obtained  merely  by  exactly  fitting 
three  constants  to  three  points,  one  each  at  high,  low,  and  intermediate 
values  of  R’ .  This  equation  is  compared  below  with  the  results  computed  by 
the  long  method  (see  table  on  page  18). 

For  practical  reasons  it  was  convenient  to  use  units  other  than  moles  and 
atmospheres.  The  sodixmi  is  given  in  hundreds  of  cubic  feet  measured  at  350°F, 
the  air  in  millions  of  cubic  feet  at  77°F,  and  the  pressure  in  pounds  per 
square  inch  gauge.  The  equation  was  used  to  draw  generalized  curves,  presented 
in  Hgure  MH-9-A1804,  which  were  sufficiently  accurate  to  be  used  for  the  final 
phases  of  the  Power  Plant  Building  designing. 

It  can  be  shorn  by  similar  reasoning  that  if  variable  water-to- sodium 
ratio  is  to  be  considered,  the  equation  for  P  becomes  of  third  degree  in  R^ 
and  Ry,  with  R^  =  air/sodiiim  and  =  water/sodium.  Since  long-method  compu¬ 
tations  had  shown  that  the  maximum  pressure  resialted  for  a  y/x  ratio  of  unity 
in  the  range  of  Rj^  shown  in  the  table,  the  water  ratio  effect  was  not  con¬ 
sidered  further. 


Comparison 

of  Eq'uation  (P_  =  3.A4- 

+  45.64E 

'  -  7.97R 

'  ■  .  - 

with.  Values  Computed  bj'- 

Long  Piethod 

Sodium 

Air 

P 

2a 

A 

lo90 

2 

0.950 

39.7 

39.6 

+  0.1 

2.23 

2 

1.115 

44*4 

44.4 

0.0 

(1)*  2.23 

2.5 

0.892 

37.8 

37.8 

0,0 

2.23 

3.0 

0.743 

33.1 

33.0 

+  0.1 

2.56 

2 

1.28 

49.1 

48,8 

+  0.3 

2.90 

2.7 

1.074 

43.3 

43.3 

0.0 

3.05 

2.5 

1.22 

47.0 

47.3 

-  0.3 

(2)*  3.50 

2.5 

1.400 

51.7 

51.7 

0.0 

3.50 

3.25 

1.077 

43*4 

43.3 

+  0.1 

3.50 

3.65 

0.959 

39.8 

39.9 

-  0.1 

3.50 

3.76 

0.931 

39.1 

39.0 

+  0.1 

3c50 

4„00 

0.875 

37.3 

37.3 

0.0 

(3)*  3.50 

6.50 

0.538 

25.7 

25.7 

0.0 

3.50 

7.00 

0.467 

22.7 

23.0 

-  0,3 

*(1),  (2) 

,  and  (3) 

were  the  points  used  to  determine  the 

equatiorio  P  is  computed  by  the  long  method,  Pg  by  the  equation, 

A  =  P  -  Pg. 

is  the  ratio  of  sodium  to  aiPo 
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